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Abstract 

We develop a sound and complete equational theory for the functional quantum pro- 



o 

O ■ 

If} , gramming language QML. The soundness and completeness of the theory are with 

respect to the previously-developed denotational semantics of QML. The complete- 
ness proof also gives rise to a normalisation algorithm following the normalisation 
by evaluation approach. The current work focuses on the pure fragment of QML 
omitting measurements. 



Key words: quantum programming, completeness, normalisation 



jj ■ 1 Introduction 



The language QML was previously introduced by the first two authors [AG04]. 
Its semantics is inspired by the denotational semantics of classical reversible 
computations. This previous work provides a semantic foundation for reason- 
ing about quantum programs by mapping them to their denotations. 

The natural next step is to develop reasoning principles on QML programs 
themselves which avoid the detour via the denotational semantics. For exam- 
ple, given the following QML definition of the Hadamard gate: 
H x = if ° x 

then (false + (—1) * true) 

else (false + true) 
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We would like to verify that H (H x) is observationally equivalent to x, using 
a derivation like: 

H (H x) = if° (if x 

then {false + (— 1) * true) 
else [false + true)) 
then [false + (— 1) * true) 
else (false + true) 

— by commuting conversion for if ° 
= if x 

then if° (false + (—1) * true) 

then (false + (—1) * true) 

else (false + tree) 
else if° (false + true) 

then (false + (—1) * true) 

else (false + tree) 

— by if 
= if x 

then (false — false + true + true) 
else (false + /a/se + true — true) 

— by simplification and normalisation 
= if ° x then true else false 

— by 7/-rule for if ° 

= x 

It is relatively easy to develop some set of sound equational principles. 
Inspired by equivalences on classical computations, one may hypothesise that 
certain equations should hold and simply verify that both sides of the equation 
have the same denotation. 

Given, however, that QML is based on a first-order functional language 
with finite types, it should be possible to also develop a complete set of equiv- 
alences that totally capture denotational equivalence. Technically, one can 
prove completeness of the equational semantics by "inverting" the denota- 
tional meaning function. The construction is subtle in parts. We present it 
first in the context of the classical sublanguage of QML, and then extend it 
to deal with quantum data and control. 

The paper is thus organised as follows. We begin with an informal review 
of QML in Section 3. In Section 4, we present the denotational semantics of 
the classical sublanguage of QML, and present a system of equations that is 
sound with respect to the denotational semantics. We then show that this set 
of equations is complete in Section 5. Section 6 repeats the development for 
the quantum constructs. Section 7 concludes. 



2 



Altenkirch, Grattage, Vizzotto, and Sabry 

2 Related work 

Peter Selinger's influential paper [Sel04] introduces a single- assignment (es- 
sentially functional) quantum programming language, which is based on the 
separation of classical control and quantum data. This language combines 
high-level classical structures with operations on quantum data, and has a 
clear mathematical semantics in the form of superoperators. Quantum data 
can be manipulated by using unitary operators or by measurement, which can 
effect the classical control flow. 

Recently, Selinger and Valiron [SV05] have presented a functional language 
based on the same classical control and quantum data paradigm. Selinger and 
Valiron's approach is in some sense complementary to ours: they use an affine 
type system (no contraction), while we use a strict system (no weakening). The 
lack of contraction is justified by the no-cloning property of quantum states. 
However, this does not apply to our approach, since we model contraction 
by sharing not by copying — this is also used in the calculus of Arrighi and 
Dowek [AD04]. 

Andre van Tonder [vT03a,vT03b] has proposed a quantum A-calculus in- 
corporating higher order programs, but no measurements. He also suggests 
an equational theory for strict (higher order) computations, but shows neither 
completeness nor normalisation. 

3 QML Syntax and Examples 

The QML terms consist of those of a first-order functional language, extended 
with quantum data and quantum control. The full language also includes 
quantum measurement, which we do not consider in this paper. The syntax 
of terms is the following: 

( Variables) x, y, ... G Vars 

(Prob. amplitudes) k, l, ... G C 

(Patterns) p,q ::= x \ (x,y) 

(Terms) t, u, e ::= x | () | (t, u) 

let p — t in u 
if° t then u else u' 
false | true | | k * t \ t + u 
The classic sublanguage consists of variables, let-expressions, unit, pairs, 
booleans, and conditionals. Quantum data is modelled using the constructs 
k * t, , and t + u. The term k * t where k is a complex number associates the 
probability amplitude k with the term t. It is convenient to have a special con- 
stant for terms with probability amplitude zero. The term t + u is a quan- 
tum superposition of t and u. Quantum superpositions are first-class values: 
when used as the first subexpression of a conditional, they turn the conditional 
into a quantum control construct. For example, if° (true+false) then t else u 
evaluates both t and u and combines their results in a quantum superposition. 
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3. 1 Examples 

To give further intuition about the semantics of QML, we consider a few 
more interesting examples. In the examples, we allow the definition and use 
of "global" function symbols. Adding such definitions to the formalism is 
possible but tedious, so we keep them at an informal meta-level. 

The following three functions correspond to simple rotations on qubits: 
qnot x = if ° x then false else true 

had x = if ° x then ((—1) * true + false) else (true + false) 
z x = if x then (i * true) else false 
The first is the quantum version of boolean negation: it behaves as usual when 
applied to classical values but it also applies to quantum data. Evaluating 
qnot (n * false + i * true) swaps the probability amplitudes associated with 
false and true. The second function represents the fundamental Hadamard 
matrix, and the third represents the phase gate. 
The function: 
cnot c x = if ° c 

then (true, qnot x) 
else (false, x) 

is the conditional-not operation, which behaves as follows: if the control 
qubit c is true it negates the second qubit x; otherwise it leaves it unchanged. 
When the control qubit is in some superposition of true and false, the result is 
a superposition of the two pairs resulting from the evaluation of each branch 
of the conditional. For example, evaluating cnot (false + true) false produces 
the entangled pair (false, false) + (true, true). 

3.2 Copying and Discarding Quantum Data 

To motivate the main aspects of the type system in the next section, we 
examine in detail the issues related to copying and discarding quantum data. 

A simple example where quantum data appears to be copied, in violation 
of the no-cloning theorem [NCOO], is: 

let x = false + true 

in (x, x) 

As the formal semantics of QML clarifies, this expression does not actually 
clone quantum data; rather it shares one copy of the quantum data. With 
this interpretation, one can freely duplicate variables bound to quantum data. 
When translated to the type system, this means that the type system imposes 
no restrictions on the use of the structural rule of contraction. 

Discarding variables bound to quantum data is however problematic. Con- 
sider the expression: 

let (x, y) = (false, false) + (true, true) 

in x 

where the quantum data bound to y is discarded. According to both the phys- 
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ical interpretations of quantum computation, and the semantics of QML, this 
corresponds to a measurement of y. Since measurement is semantically quite 
complicated to deal with, we insist that it should be represented explicitly. 
The language we consider in this paper lacks the explicit constructs for mea- 
surement so we reject the expression above. This means that the structural 
rule of weakening is never allowed in situations where information may be lost. 

4 The Classical Sublanguage 

By the classical sublanguage, we mean the subset of terms excluding quantum 
superpositions and hence quantum control. 

4-1 Type System 

The main role of the type system is to control the use of variables. The typing 
rules of QML are based on strict linear logic, where contractions are implicit 
and weakenings are not allowed when they correspond to information loss. As 
explained in the previous section, weakenings correspond to measurements, 
which are not supported in the subset of the language discussed in this paper. 

We use a, r, p to vary over QML types which are given by the following 
grammar: 

o = Q\ I 0,2 I o (g> r 
As apparent from the grammar, QML types are first-order and finite: there 
are no higher-order types and no recursive types. The only types we can 
represent are the types of collections of qubits. 

Typing contexts (T, A) are given by: 
r = • 1 r, x : a 

where • stands for the empty context, but is omitted if the context is non- 
empty. For simplicity we assume that every variable appears at most once. 
Contexts correspond to functions from a finite set of variables to types. We 
introduce the operator ®, mapping pairs of contexts to contexts: 

(T, x : a) <g) (A, x : a) = (T ® A), x : a 

(r, x : a) <g> A = (r <g> A), x : a if x (£ dom (A) 

• <8> A = A 

This operation is partial: it is only well-defined if the two contexts do not 
assign different types to the same variable. Whenever we use this operator we 
implicitly assume that it is well-defined. 

Figure 1 presents the rules for deriving valid typing judgements r h t : a. 
The only variables that may be dropped from the context are the ones of 
type Q\ which, by definition, carry no information. Otherwise the type system 
forces every variable in the context to be used (perhaps more than once if it 
is shared). 
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var 



rh£:<7 A, x : a h u : r 

let 



x:crhx:a r® Ah let i = t in ii:t 

r h i : a Ah«:r 



unit 



h : 2i r® Ah (i,u) :a(g)T 

r h i : a <S) t A, x : cr,y : t h u : p 

<g 

r®Ah let (x,y) = t in u : p 



-intro 



-elim 



f-intro t-intro 

h false : Q2 • h true : Q2 



Thc:Q 2 Aht,u:a T,x : Q x \- t : a 

if wk-umt 



r (8> A h if c then t else u : cr T h t : 



Fig. 1. Typing classical terms 
4-2 The Category of Typed Terms 

The set of typed terms can be organised in an elegant categorical structure, 
which facilitates the proofs later. The objects of the category are contexts; 
the homset between the objects T and A, denoted TrnTA, consists of all the 
terms t such that r h t : |A| where |A| views the context A as a type. This 
latter map is naturally defined as follows: 

I • I = Qi 

\T,x : a\ = \T\ ® a 
For each context T, the identity lr € Tmrr is defined as follows: 

l. = () 

lr,x:o- = (lr, %) 

To express composition, we first define: 
let* • = u in t = t 

let* r, x : a = u in t = let (x r , x) = u in let* r = x r in t 
Given d 6 Tm A T and e G TrnT 0, the composition e o d G Tm A is given 
by the term let* r = d in e. 

4-3 Semantics 

The intention is to interpret every type a and every context T as finite sets [a] 
and [r], and then interpret a judgement r h t : a as a function |F h t : cr] G 

[r] - H- 
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In the classical case, the type Q 2 is simply the type of booleans; the types 
are interpreted as follows: 

[Qi] = {0} 
[Qa] = {0, 1} 
fa S r] = M x [r] 

We use the abbreviation [r] for [|r|J. 

The meaning function is defined in Figure 2 by induction over the structure 
of type derivations. It uses the following auxiliary maps: 

• id : S — > S defined by id (a) = a 

• id* : S — > [Qi] x S and its inverse id* defined by id* (a) = (0, a) and 
id*(0, a) = a 

• For a E S, the family of constant functions const a : [<2i] — > S defined by 
(const a)(0) = a. 

• 5 : S — > (S, S) defined by 5(a) = (a, a) 

• swap : S x T — > T x S defined by swap(a,b) = (b,a). We will usually 
implicitly use swap to avoid cluttering the figures with maps which just 
re-shuffle values. 

• For any two functions / G Si — > T\ and g G S2 ~^ T 2 , the function (/ x g) : 
(Si x S2) — > (Tx x T 2 ) is defined as usual: 

(/ x g)(a,b) = (f a,g b) 

• 5r,A ■ [r <S> A] — > [r] x [A]. This map is defined by induction on the 
definition of T <S> A as follows: 



Sr',A' x S if T = T', x : a and A = A', x : a 
Sr',A x id if F — T', x : a and 2 G" dom (A) 
id* if T = • 



Intuitively, the map S-p,A takes an incoming environment for an expression, 
creates shared copies of the appropriate values, and rearranges them (the 
shuffling is implicit and not shown in the above definition) into two envi- 
ronments that are then passed to the subexpressions. 

For any two functions /, g G S — > T, we define the conditional f\g G 
([Q2I x S) -> T as follows: 

(f\g) (l,a) = fa 
(f\g) (0,a)=ga 
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ir« u a . 
r H • 


Will 


= const 


[[• r raise . 


D H 


= const 


[• h true : 


W^J 


= const 1 


[i : ff h i 


■ ^-H 
• 


.'j 

— ?u* 


[T"P a A ^ „ . 

[[i cx) za i ieu x — i in u 


• T J 


= 9 ° (f x id) o 5r,A 






wnere / = |]1 r I : oj 






9 — H^*) x . a r u . rj 


|[l (x) za r ^c, MJ . cr c 




— [J x 9) ° °T,A 






wnere j = |[1 r r : uj 






9 — h u . rj] 


H i (xi za r lex \x, y ) — t m ii 


• PI 


= 9 ° (f x id) o Sr,A 






wnere / — |[i r i . <r t>9 tj 






5 = [A, x : o-,y : t h u : p} 


ii x " r 11 c tnen r else u 


• 0J 


— U?I"J ° U x %a ) ° »r,A 






wnere J - [l r c . Wi2jJ 






5 = [A h t : a] 






/i = [A h u : er] 


[Tht 




= foid* 






where / = [T, x : Qi h i : cr] 



Fig. 2. Meaning of classical derivations 



^.^ Equational Theory 

We present the equational theory for the classical sublanguage and then show 
its soundness and completeness. The equations refer to a set of syntactic 
values defined as follows: 

val G Val c ::= x | () | false \ true \ (val\, vah) 

Definition 4.1 The classical equations are grouped in four categories. 

• let-equation 

let p = val in u = u [val / p] 

• /3-equations 

let (x, y) = (t, u) in e = let x = t in let y = u in e 

if ° false then t else u = u 

if° true then £ else u = t 

• ^-equations 
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t — if t:Q 



let x = t in x 

let (x, y) = tin (x, y) 

if ° t then true else false 



t 
t 
t 



• Commuting conversions 

let p = t in let q = u in e 



let q = u in let p = t in e 
if t 

then let p = Uq in e 
else let p = u\ in e 



let p 



if ° i 
then «o 
else ui 



in e 



We write rht = u:crifri-t,u:<7 and the equation t = u is derivable 
at the type a. 

Lemma 4.2 (Soundness) TTie equational theory is sound: iJY \- t = u : a 
then the functions [r h t : <r] and |r h u : a] are extensionally equal. 

5 Completeness of the Classical Theory 

The equational theory is complete in a strong technical sense: as we prove 
in the remainder of the section, any equivalence implied by the semantics is 
derivable in the theory. The proof technique is based on current work by the 
first author with Tarmo Uustalu [AU04]. The proof we present extends and 
simplifies the method presented in that work. 

5.1 Proof Technique 

The ultimate goal is to prove the following statement. 

Proposition 5.1 (Completeness) If \T h t : aj and [r h u : er] are exten- 
sionally equal, then we can derive V h t = u : a. 

In order to prove this statement, we define a function q£ which inverts 
evaluation by producing a canonical syntactical representative. In fact, we 
define the function such that it maps a denotation [T h t : er] to the normal 
form of t. 

Definition 5.2 The normal form of t is given by nf£(t) = Qr([r h i : er]). 

The normal form is well-defined: given an equation r h t = u : a, we know 
by soundness that [r h t : a] is extensionally equal [r h u : a] and hence we 
get that nfp(t) = nf^(it). If we can now prove that the syntactic theory can 
prove that every term is equal to its normal form, then we can prove the 
main completeness result. Indeed given the following lemma, we can prove 
completeness. 

Lemma 5.3 (Inversion) The equation V h nf£(t) = t : er is derivable. 
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Proof of Proposition 5.1 (Completeness) We have: 

r h t = c/f |T h t : cr] : o by inversion 
T h gf[r h t : <tJ = gf[r h n : o"J : <T by assumption 
r h c/f [r h u : cr] = u : a by inversion 

□ 

To summarise we can establish completeness by defining a function that 
inverts evaluation and that satisfies Inversion Lemma 5.3. 

5.2 Adequacy 

We begin by defining a family of functions q a ("quote") which invert the 
evaluation of closed terms and prove a special case of the inversion lemma for 
closed terms, called adequacy. These functions and the adequacy result are 
then used in the next section to invert the evaluation of open terms and prove 
the general inversion lemma. 

Definition 5.4 The syntactic representations of denotations is given by: 

q a G H -> Val c a 

defined by induction over a: 

q Ql = () 
c/ 2a = false 
q Q2 1 = true 
q °®T ( a , b) = (q a a, q T b) 

The version of the inversion lemma for closed terms is called adequacy. It 
guarantees that the equational theory is rich enough to equate every closed 
term with its final observable value. 

Lemma 5.5 (Adequacy) The equation h q a (\ h t : cr] 0) = t : a is deriv- 
able. 

Proof sketch. During the proof of such a statement we encounter open terms 
that must be closed before they are "quoted." So in fact the statement to prove 
by induction is the following: 

If g G [r] then h <f([r h t : aj g) = let* r = q T (g) in t : a 

□ 

5.3 Inverting Evaluation 

As explained earlier, the main ingredient of the proof of completeness is the 
function c/f which inverts evaluation. To understand the basic idea of how the 
inverse of evaluation is defined, consider the following example. Let T be the 
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environment x : (Q 2 <S> Q 2 ), V '■ Q2 and let / G [r] — > [22]- To find a syntactic 
term corresponding to /, we proceed as follows: 

• flatten all the products by introducing intermediate names; this produces 
an updated environment V — x± : Q 2 ,x 2 '■ Q.2,V '■ Q2, and an updated 
semantic function /' such that: 

/'(((((), x 1 ),x 2 ),y) = f((0,(x 1 ,x 2 )),y) 

• enumerate all possible values for the variables, and apply /' to each enumer- 
ation to produce a result in the set 10,2}- For example, it could be the case 
that / (((), (1, 1)), 1) = 0. The result of each enumeration can be inverted 
to a syntactic term using q a from Definition 5.4. 

• Put things together using nested conditions representing all the possible 
values for the input variables. In the example we are considering, we get: 

let (xl, x2) = x 
in if xl 

then if x2 

then if y then false 
else 

else ... 
else ... 

The idea is formalised in the following definition. 
Definition 5.6 The function 

for inverting evaluation is defined by analysing the context: 

<£(/) = <f (f (0)) 

«?«*:&(/) = It (h) where h{g) = f(g, 0) 
Qr, X :Q 2 (f) = ( if° x then (h) else {ho)) 

where h^g) = f(g,i) for i e {0, 1} 
2rWi®r 2 )(/) = ( let {xi,x 2 ) = x in q^ Xi:TuX2:T2 (h) 

where h(g,x 1 ,x 2 ) = f(g,(xi,x 2 )) 

The base case is straightforward: the evaluation produces a closed value 
which can be inverted using the "quote" function of Definition 5.4. If the 
context includes a variable x of type Qi, then we supply the only possible value 
for that variable (0), and inductively construct the term with the variable x 
bound to (). The result is of the correct type because we can add or drop 
bindings of variables of type Q\ to the environment. If the context includes a 
variable x of type Q 2 , then we supply the two possible values for that variable 
and 1. A conditional is then used to select the correct branch depending on 
the actual value of x. Finally, if the context includes a variable of type t\ <8> t 2 
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z-intro 



r h t : a 



r h t,u : a 



• hO:(T 



prob 



Vrt + u-.cr 



sup 



rh K*t : a 



Fig. 3. Typing quantum data (I) 

then we simply flatten the product and proceed inductively. The function q£ 
does indeed satisfy the inversion lemma. 

6 Quantum Data and Control 

We develop the typing rules and semantics of the quantum fragment of QML 
in two stages. First we extend the judgements T h t : a and the semantics of 
Section 4 to handle quantum data in a straightforward manner. This simple 
treatment is only however an intermediate step in the development as it admits 
quantum programs that are not realisable on a quantum computer. We then 
refine both the type system and the semantics to identify exactly the realisable 
quantum programs. 

6. 1 The Category Vec 

As a first approximation to a type system for QML programs, we consider the 
type system of Figure 1 extended with the rules in Figure 3. 

Unlike the classical judgement Y h t : a is not interpreted as 

a function in [F] — > [er]. Rather, because we now have superpositions of 
terms with complex probability amplitudes, we interpret such judgements as 
functions in [T] — > [erj^ where [cr]^ represents the complex vectors over the 
base set [cr]. In other words, [crj^ is defined to be [<r] — > C which is sometimes 
denoted V [cr]. We call the structure described above the category Vec. 

Naturally this change requires that we revisit the semantics of the classical 
terms given in Figure 2 so that each denotation returns a complex vector. For 
example, we should have: 



Instead of mapping the value representing the empty context to the deno- 
tation of false, we now return a vector v which associates the denotation of 
false with probability amplitude 1 and the denotation of true with probability 
amplitude 0. 

This change can be done systematically by noticing that it corresponds to 
a monad whose unit and lift operation are defined below: 



More precisely every value that is returned in Figure 2 is explicitly tagged 
with the monadic return and when two functions are composed in Figure 2 
using fog, the composition is replaced by /* o g. 



[• h false : Q 2 j 



Q 



const v where v = 1 and v 1 = 



return a (b) = 1 if a = b and otherwise 
f*(v) = Za.(v a) * (/ a) 
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[r h k * t : aj Q 



const v where Va G [cfl.u a = 
g where g a = k* (fa) 



IT h t + u : a] Q 




Fig. 4. Meaning function for quantum data 



The meaning of the new constructs for quantum data is given in Figure 4. 



6.2 Orthogonality 

The type system presented so far does indeed correctly track the uses of vari- 
ables and prevents variables from being weakened; yet the situation is more 
subtle. It turns out that the type system accepts terms which implicitly per- 
form measurements and as a consequence accepts programs which are not 
realisable as quantum computations. 

Consider the expression if° x then true else true: this expression appears, 
syntactically at least, to use x. However given the semantics of if , which 
returns a superposition of the branches, the expression happens to return 
true without really using any information about x. In order to maintain the 
invariant that all measurements are explicit, the type system should reject the 
above expression as well. 

More precisely, the expression if° x then t else u should only be accepted 
if t and u are orthogonal quantum values (t _L u). This notion intuitively 
ensures that the conditional operator does not implicitly discard any infor- 
mation about x during the evaluation. Because of a similar concern, the two 
branches of a superposition should also be orthogonal. 

The typing rules for conditionals and superpositions are modified as in 
Figure 5. This modification also achieves that programs are normalised, i.e., 
the sum of the probabilities of a superposition add up to 1. 

In Figure 6 we define the inner product of terms, which to any pair of 
terms r h t, u : a assigns (t\u) G C U {?}. This is used to define orthogonality: 
t JL u holds if (t\u) = 0. 

The judgement h° is not automatically closed under the equality judge- 
ment, hence we add the rule (subst). Our philosophy is that we allow equiva- 
lent representations of QML programs which do not satisfy the orthogonality 
criteria locally, as long as the program as a whole is equivalent to one which 
does satisfy the criteria. 
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r h° c : Q 2 Ah°t,u:a t±u 

if 

T<g)AI- if° c then t else u : a 

Th°t,u:a t-Lu |A| 2 + \k\ 2 = 1 

si 

F\-°X*t + K*u: a 

r h° t:a Tht = u:a 

subst 

T h° u : a 



Fig. 5. Typing quantum data (II) 



(tit) 


= 1 




(A * t + A' * t! | u) 


= X* * (t\u) + X H * (t'\u) 


(false true) 


= 




(t K * U + k' * u') 


= K * (t\u) + K 1 * (t\vf) 


(true false) 


= 














(X*t\u) 


= X*(t\u) 


(If (true) 


= = 


(true| If) 


(t\X*u) 


= X(t\u) 


(0*1 false) 


= = 


(false) If) 


(t + t'\u) 


= (t\u) + (t'\u) 


(oV) 


= = 




(t\u + It') 


= (t\u) + (t\u>) 


<(t,f) | («,«')> 


= (t\u 


} * (t'W) 


(t\u) 


= ? otherwise 



Fig. 6. Inner products and orthogonality 



6.3 The Category Q° 

The restriction of the set of typable terms requires a similar semantic restric- 
tion. All we need to do is to restrict the morphisms in the category of complex 
vectors to satisfy the following two conditions: 

• Linearity: If / G V A — > V B, a G C, and v, v±, v 2 G V A, then f(vi + v 2 ) = 
f(Vi) + f(v 2 ) and f(av) = a(f v). 

• Isometry: If/eVi-*VB and Vi,v 2 G V A, then (vi\v 2 ) = (f Vi\f v 2 ). 
(In other words, / preserves inner products of vectors.) 

Two morphisms f,g G A — >• B are orthogonal if for all vector v G V A, we 
have (/ v\g v) = 0. We call the resulting category, the category Q° of strict 
quantum computations. The homset of morphisms in [r] — > [<r]p satisfying 
the above conditions is called Q° [r] 

The meaning function is given as before but with the maps interpreted in 
the category Q°, i.e., the meaning of a derivation r h t : a is a morphism 
[r h t : cr] Q G Q° [r] [<r] Q . The requirement for orthogonality in the type 
system is reflected semantically: for isometries f,g, we have that f\g is an 
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isometry, if / and g are orthogonal. 

6.4 Quantum Equational Theory 

The equational theory for the quantum language inherits all the equations 
for the classical case. This can be informally verified by noting that the 
meaning function in the case of the quantum language is essentially identical 
to the classical case. Formally, the proof technique explained in Section 4 
applies equally well to the quantum case and yields the same equations for 
the classical core plus additional equations to deal with quantum data. 

Definition 6.1 The quantum equations are: 
(if) 

if (A * to + k * ti) then uq else m 
= A * (if° to then uq else ui) + K * (if h then no else U\) 

(superpositions) 

t + u = u + t 

t + 1? = t 

t + (u + v) = (t + u) + v 

A * (t + u) = \ * t + \ * u 

A * t + k* t = (A + k) * t 

0*t = 7f 

Lemma 6.2 (Soundness) The equational theory is sound: i/Th t = u : a 
then the isometries [r h t : a] ® and [r h u : cr] ^ are extensionally equal. 

The additional equations are used to prove equality between different quan- 
tum values. Semantically, two quantum values are the same if they denote the 
same vector, which is the case if the sum of the paths to each classical value 
is the same. For example, to find a simplified quantum value equivalent to: 

{false + true) + [false + (— 1) * true) 
we first normalise to: 

(1 / y/2) * ((1 / \/2) * false + (1 / V2) * true) + 

(1 / V2) * ((1 /y/2)* false + (-1 / V%) * true) 
This term has two paths to false; along each of them the product of the 
amplitudes is (1 / V2) * (1 / V2) which is 1 / 2. The sum of all the paths 
to false is 1, and the sum of all the paths to true is 0. In other words, the 
entire term is equivalent to simply false. The above calculation proves that 
the Hadamard operation is self-inverse, as discussed in the introduction. 

6.5 Quoting quantum values 

We will now adapt the techniques developed in section 4 to the quantum case. 
A classical value v £ Val c <7 is simply a term representing an element in [a]. 
A quantum value represents a vector in V [er]p, hence we have to close values 
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under superpositions. We define Val^ a C Tm a inductively as a subset of 
closed terms of type a: 

v G Val c a 

• 

val v G Val^ a 

• G Val Q a 

v, w G Val Q a 

• 

v + w G Val Q a 
v G Val Q a 

• 

« * v G Val Q cr 

We write Valuer for isometric quantum values which satisfy the restrictions 
introduced in Figure 5. 

We have already seen that there is a monadic structure on V A = A — > C. 
Correspondingly, we have a Kleisli structure on Val^; val G Val c <7 — > Valuer 
is the return and bind is defined as given v G Val Q cx and / G Val c a -> Val Q r, 
we define v f G Val^ r by induction over v: 

(val x) f = f x 

^/ = 

f + w »= / = (v ^= f) + (w >= /) 

ft * V ^=f = K*(v^^f) 

Lemma o.o (Val c ,Val Q ,wZ, (»=)) «5 a Kleisli structure, i.e. it satisfies the 
following equations: 

(i) wa/ a; f = f x 

(ii) v 3= Xx.val x = v 

(hi) v »= Ax.(/ z) »= g = (v »= /) »= # 

Proof. Case (i) follows from the definition. Cases (ii) and (iii) can be shown 
by induction over the structure of v. □ 

While the classical definition of q a (def. 5.4) was completely straightfor- 
ward, its quantum counterpart is a bit more subtle, in particular the in the 
case of tensor products. As a special case consider g Q2 ®2 2 ; given an element 

IT e\Q 2 ® Q 2 J Q = [Q 2 ] x lQ 2 j -> C 

we have to construct a value q Q2 ® Q2 if g Val^ 0,2® 0,2- This can be done by 
calculating the probabilities that the first qubit is i, fst ~v*i G M + , given by 

fst If i = v/|l?(i,0)| 2 + \l?(i, 1)| 2 

creating the first level of the value as a tree, and then for the second level 
normalising the amplitudes wrt. the probabilities of the previous level, see 
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"(0,0) / 


\ "(0,1) 


v(l,o) y 


\ "(1,1) 


fst v / 


\fst v 


fst V l ./ 


\f«t v 1 


(0,0) 


(0,1) 


(1,0) 


(1,1) 



Fig. 7. Value tree for Q 2 (8) Q2 

figure 7 for the corresponding tree. We write [<r] p = [<r] — > M + for the set 
of probability distributions, obviously we have [<t] p C [a] Q . We observe that 
fst if G [crj p . Generalising the idea given above we arrive at the following 
definition of quote: 

Definition 6.4 The syntactic representations of denotations is given by 

<f G H Q -> Val Q a 

defined by induction over a: 

qQi if = (if 0) * () 

qQ.2 -iff = (-if 1^ % ^rue _|_ (jf q-) * f a i se 

^\x G |cr].(l/(fstlf)x) * q T (Xy.lf(x,y)) 
^=\y.val (x, y) 

where: 

fst G [a ® r] Q -> [a] p 
fst 7? x = ^/Ey.|lf (x, y)\ 2 

l/_ e M p _ M p 

l/lfx = Ax. if p^eO then else 1 / (p 1) 

To show adequacy we have to establish a number of properties of q a : we 
have to show that it is linear and isometric and that it preserves tensor prod- 
ucts. This is summarised in the following proposition: 

Proposition 6.5 

(i) q a (k * if) ee k * (q a if) 

(ii) q a (if + Uf) = (q a if) + (q a if) 

(iii) (if I vf) = (q a lf\q a vf) 

(iv) q a ® T (!f<8>vt) = (q a if, q T vf) 

The proof of the above proposition again isn't completely straightforward, 
e.g. linearity cannot just be proven by induction over a. It is essential that 
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we first establish some properties of renormalising a vector wrt. a probability 
distribution. We define the product of a probability distribution p G [cr] p and 
a vector if G [foi] Q as: 



p * v G hr 



IQ 



p * if = \x G [cr].(px) * (if x) 

It is not hard to see that an analogous operation can be defined on values, 
given v G Val^ a and p G [<r] p as above, we define: 

p * v G Val Q a 

p * v — v >■= \x G [o"J.(px) * (valx) 
The key property we establish is 

Lemma 6.6 Given p G [cr] p and if G [<r]p 

p * if) = q a (p * if) 

which can be verified by induction over a and observing that while 1/— 
isn't a proper inverse, it nevertheless satisfies the following property 

l/(p + q) *(p + q) = (l/p)*p 

Using the fact that q a is isometric we can show that it produces values 
satisfying the orthogonality constraints: 

Proposition 6.7 Given v G [cr] <3 

h° q a v : a 



6.6 Adequacy 

We define a syntactic counterpart to: 

5r,A e Q° [r®A]([r]Q®[A] ( 



as: 



by: 



5 r ,A GTm(r® A) (|r| ® |A|) 



5 T A = < 



let ((/, d) = <5r',A'i n ((g, x ),(d, x)) if Y — V, x : a 
and A = A', x : a 
let d) = <5r',Ain ((^5 d) if T = T' , x : a 
and x G" dom A 

1 A if r = • 

To establish that q a commutes with the context operations we have to 
show that contraction corresponds to S G Q° [er] (|cr] Q <S> |o r ] Q )- 
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Lemma 6.8 Given v E [er] ® we have 

let x = q a v in (x, x) = t> 

Proof. By induction on er. □ 

Exploiting this property we can show that the context operations commute 
with quote: 

Lemma 6.9 Given if E [r <g> AJ Q 

g |r|»|A| {6rA ^) = s rA q^ if 

Theorem 6.10 IfT\-t:a and g E [T] Q then 

h <f([r h t : cr] Q #) = Zef T = q r g in t : <r. 

Proof. By induction over the derivation of V h £ : er, as an example consider 
the case for let: 

q p ([r <g> A h let x = t in u : p] Q ) 

= {definition of [. . .J Q } 

qP o (|t]Q <g> id) o 5r )A ) 

= {induction hypothesis for u and £} 

M o (£ o g r ® g A ) o 5 rA ) 

= { lemma 6.9 }} 

u o (£ ® id) o <5 r> A o gl r ® A l 



(let x = £ in u) o g 



|r®A| 



The other cases use the same style of reasoning to deal with the structural 
properties and exploit proposition 6.5. Note that the case for if° can be 
reduced to linearity. □ 

Corollary 6.11 (Adequacy) J/h £ : er then hf([ h t : oj Q ) = t : er 



6. 7 Completeness and normalisation 

The development here follows closely the one in the classical case as presented 
in Section 5.3. 

Definition 6.12 The function: 

e Q° [r] [a] Q ->Tmrer 
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for inverting evaluation is defined by analysing the context: 

q%(f) = (f (/ (return 0)) 

#,*:&(/) = ^Fi:Qi ° (?r) ° $ r,x: Sl 
Qv,x:Q 2 {f) = 0F,U 2 ° (?f X 9?) ° $ r,x:fi 2 

The auxiliary isomorphisms are defined as follows: 
0r> : Q! <E Tm (r, x : Q\) a — > Tm T a 
rfr^-.Q^ = let x = () in t 
<M = * 



0r,x:Q 2 G Tm(r,x : Q 2 a) -> {(t , *i) e (TmlV) 2 | t -L ti} 
0x:Q 2 ^ = (let ^ = false in t, let x = true in t) 
0P^,.Q 2 (t,M) = if x then t else u 



<fa,xvn&n e Tm (r, x : n <g> r 2 ) p -> Tm (r, Xi : n, x 2 : r 2 ) 

<t>v,x:r 1 ®T2 t = let x = (xi,x 2 ) in t 

0r,x:Ti^(*) = let {xi,x 2 ) = x in t 
The semantic map corresponding to each is written $. 

For the inversion proof we only need the provability of one side of the 
isomorphisms which follows from the //-equalities. 

Lemma 6.13 The following family of equalities is derivable 

^ l ((j) T t)=t 

Definition 6.14 The normal form of t is given by nf£(t) = ([r h t : crj^). 
Lemma 6.15 (Inversion) The equation V h nf£(t) = t is derivable. 

Proof. By induction over the definition of q£- In the case of T = • the result 
follows from adequacy, Corollary 6.11. In all the other cases we exploit Lemma 
6.13. □ 

Since all our definitions are effective nf indeed gives rise to a normalisation 
algorithm. As a consequence, our equational theory is decidable, modulo 
deciding equalities of the complex number terms which occur in our programs. 
We also note that as in the classical case, our theory is complete: 

Proposition 6.16 (Completeness) If |T h t : a] Q and [r h u : a} Q are ex- 

tensionally equal, then we can derive V h t = u : a. 

20 



Altenkirch, Grattage, Vizzotto, and Sabry 

7 Conclusions and Further Work 

We have developed a sound and complete equational theory for a functional 
quantum programming language, while at the same time providing a nor- 
malisation algorithm. The construction is a modular extension of a classical 
theory, indeed the quantum theory inherits not just all the equations and term 
formers, it is also possible to generalise our proof technique to the quantum 
case. The quantum theory introduces additional constructs corresponding to 
superpositions and equations relating them. 

The obvious next step is to generalise this approach to the full language 
QML including measurements. The equational theory is already a challenge, 
since a measurement can have non-local effects on shared data. Semanti- 
cally, we will be using superoperators to model programs with measurements. 
Clearly, we have to extend our quote operator to work on density matrices. 

Another interesting direction, would be to consider higher order quantum 
programs and develop a complete equational theory and normalisation al- 
gorithm for this calculus. A likely semantic domain is given by presheaves, 
here the tensor product can be modelled using Day's construction, which is 
automatically closed, i.e., provides an interpretation for higher types. 
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